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Reconciling optical and X-ray mass estimates: the case of 
the elhptical galaxy NGC3379 
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ABSTRACT 

NGC3379 is a well studied nearby elliptical for which optical investigations have 
claimed a little dark matter content, or even no dark matter. Recently, its total mass 
profile M(r) has been derived by exploiting Chandra observations of its extended and 
X-ray emitting interstellar medium, based on the hypothesis of hydrostatic equilibrium 
for the hot gas. The resulting total mass within the effective radius Rc has been claimed 
to be a few times larger than found by optical studies. Here we show that part of 
the discrepancy can be due to an underestimate of the optically derived mass, and 
the remaining discrepancy of a factor of '--^ 2 can be explained by deviations from 
hydrostatic equilibrium of the hot gas. By using hydrodynamical simulations tailored 
to reproduce the observed hot gas properties of NGC3379, and by assuming as input 
for the simulations the total mass profile derived optically, we show that i) the hot 
gas at the present time has X-ray properties consistent with those observed only if 
it is outflowing over most of the galactic body, and ii) an overestimate of M of the 
same size found in the recent X-ray analysis is recovered when assuming hydrostatic 
equilibrium. We also show that the hot gas is outflowing even for a dark matter 
fraction within Rc as large as derived with the standard X-ray procedure based on the 
hydrostatic equilibrium assumption, which shows the unapplicability of the method for 
this galaxy. Finally, we find that the whole range of dark mass amount and distribution 
allowed for by optical studies is compatible with a hot gas flow with the observed X-ray 
properties. 

Key words: galaxies: elliptical and lenticular, CD - galaxies: individual: NGC3379 
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ISM 



1 INTRODUCTION 

The presence of an extended and X-ray emitting interstellar 
medium (ISM) in early type galaxies has been often used as 
a tool to determine their dynamical properties, in particular 
their mass profile and total mass out to large radii (~ 10 op- 
tical effective radii -Re; e.g., Fabbiano 1989). Most recently, 
this kind of analysis has been performed for a sample of 53 
early type galaxies by exploiting Chandra data (Fukazawa 
et al. 2006, hereafter F06). In the commonly adopted ap- 
proach, a spherical mass distribution is assumed and, in the 
hypothesis of hydrostatic equilibrium, the total mass profile 
is recovered via the well known formula (Fabricant, Lecar & 
Gorenstein 1980): 



G/irrip \ dlogr dlogr J 



(1) 



where T and p are the radial temperature and density pro- 
files, inferred from deprojection techniques (e.g.. White, 



Jones & Forman 1997), 



is the proton mass, p is the 



mean molecular weight and k is the Boltzmann constant. 
Using eq. (1), several authors measured the total mass of 
X-ray bright, nearby early type galaxies, finding in general 
larger masses than obtained with optical studies (e.g., Ciotti 
& Pellegrini 2004 and references therein; hereafter CP04) , or 
sometimes lower than these, in the central galactic regions 
(Tsai 1994). In fact mass profiles can be derived also from 
the observed stellar velocity dispersion profiles ap{R); these 
extend out to ~ 1 — 2_Rc in few best cases (e.g., Berlin et 
al. 1994, Kronawitter et al. 2000). For example, by solving 
the Jeans equations for assumed visible and dark density 
profiles and a specified radial trend of anisotropy in the 
stellar velocity dispersion tensor, one obtains o-p{R) to be 
compared with the observed one. This method suffers from 
a "dark matter-anisotropy degeneracy", because radial or- 
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bital anisotropy can mimic the presence of a massive dark 
matter halo (e.g., Tonry 1983). This technique has also been 
used considering as tracers planetary nebulae or globular 
clusters (Romanowsky et al. 2003, Teodorescu et al. 2005; 
Pierce ct al. 2006) , giving a variety of results. However, opti- 
cal studies converge towards an amount of dark mass that is 
less than the optical mass within Re [or even much lower, see 
Cappellari et al. (2006) for the most recent analysis based 
on SAURON data], and possibly larger outside. The use of 
the hot ISM to derive dynamical quantities using cq. (1) is 
then appealing both for its capability to reach large radii 
and for being exempt from the orbital anisotropy degener- 
acy problem. 

In a recent work, however, the accuracy of the X-ray 
based estimates of the total mass with cq. (1) [hereafter 
M^t(r)] has been questioned (CP04). Instead of eq. (1), a 
more general expression for the momentum equation of the 
hot ISM is (e.g., Ciotti et al. 1991): 

k dpT 



(2) 



Urupp dr 

where g = GM{r)/r'^ describes the total gravitational field 
and V is the kinetic term that includes all the non hydro- 
static contributions, such as the time dependence of the ve- 
locity, the advective term and the presence of mass sources. 
CP04 investigated the effect that deviations from hydro- 
static equilibrium have on M^((r) derived via eq. (1), and 
showed that 

2 

(3) 



They found that less massive distributions than true are 
deduced for the central regions of galaxies hosting inflowing 
gas, and conversely that the true mass is overestimated for 
outflowing regions where the velocity increases with radius, 
as in galactic winds. 

Currently, quite accurate T and p profiles are being 
derived from Chandra and XMM — Newton observations 
(e.g.. Sun et al. 2003; F06). F06, in particular, investigated 
the dynamical structure of 53 elliptical galaxies, both X-ray 
bright and X-ray faint; by exploiting the high angular reso- 
lution of Chandra, they obtained a temperature and density 
profile from 0.1 to 10 Re. For the well studied nearby ellip- 
tical NGC3379, they derived an of ^ 7 times the value 
obtained with optical methods, for the same galactic region. 
Here we demonstrate how deviations from hydrostatic equi- 
librium can be responsible for the discrepancy between M^^t 
and the mass estimated from optical studies, for this galaxy. 
By using hydrodynamical simulations tailored to reproduce 
the hot gas properties observed by Chandra for NGC3379, 
it is shown that an overestimate is obtained from the use 
of eq. (1), since the hot ISM is largely outflowing. In Sec- 
tion 2 we briefly describe the previous mass estimates for 
NGC3379 coming from optical studies and from the recent 
X-ray analysis; in Section 3 we present the galaxy model 
used for the hydrodynamical simulations and the results of 
the simulations; in Section 4 the conclusions are presented. 



Vaucouleurs (1948) law over a span of 10 magnitudes (Ca- 
paccioli et al. 1990; see Table 1 for the main galaxy proper- 
ties). Its optical shape is roundish (its morphological type is 
El) and its optical colors are typical of an old population; it 
shows no sign of current or past interactions (Schweizer & 
Seitzer 1992). The total mass profile of NGC3379 has been 
studied several times using optical data (hereafter, all quan- 
tities involving distance have been rescaled for the distance 
in Table 1; M indicates the total mass, M« the stellar mass 
and Mh the dark halo mass). 



2.1 Results of optical studies 

Even when performed with different approaches, optical 
studies invariably show very little dark matter within Re. 
Their results arc usually expressed in terms of the B-band 
total mass-to-Ught ratio M/Lb. For NGC3379 the contri- 
bution of M^/Lb to the total M/Lb has been estimated 
to range from 4.2 to 8.8 from stellar population synthesis 
models (Gerhard et al. 2001). 

The dynamical modeling by Saglia, Bcrtin & Stiavelli 
(1992), based on kinematical data extending out to ^ 1.7Re, 
gives a best fit M/Lb = 7.8 for a one-component mod- 
cling, for the whole galaxy. In their phase-space based two- 
component modeling, they found AI,/Lb = 7.5, and a dark- 
to- luminous mass ratio Mh/Mt = Q.67t.o'\o within Re, and 
Mh/M^ = 1.4to;7 for the whole galaxy. 

Kronawitter et al. (2000) 's dynamical modeling, based 
on stellar kinematic information including line profile, and 
extending from 0.005 to 2.75i?c, shows an outward increase 
of M/Lb of a factor of 1.21q'2 > with respect to a central value 
of M/Lb = 5.6^0^37 (though the modeling is also consistent 
with a flat M/Lb over the range of the kinematic data). At 
7?othey give^ M/Lb = 6.3, which implies Mh/M, = 0.13 
when for M^/Lb the central M/Lb is adopted. At their 
outermost datapoint, M^/Mt, = 0.2~^q 2- 

More recent modeling of long-slit data extending to 
80" found again that dark matter is not required out to 
~ 2.1_Rc (Samurovic & Danziger 2005); the data can be fit- 
ted with a constant M/Lb = 6 — 11. SAURON observations 
of the central region of NGC3379 give integral-field stel- 
lar kinematics out to 28" (Cappellari et al. 2006). Coupled 
with 2-integral Jeans and 3-integral Schwarzschild dynami- 
cal models, the data arc consistent with a constant mass to 
light ratio in the I band; from their M^/Lj for the stellar 
population, Mh/M^ ~ 0.09 is recovered. 

The mass profile has been studied also at larger ra- 
dial distances using planetary nebulae as tracers. CiarduUo, 
Jacoby, & Dejonghe (1993) presented a radial velocity sur- 
vey of planetary nebulae with galactocentric radii <4:.5Re. 
At large radii, the stellar dispersion profile joins with that 
of the planetary nebulae. The velocity dispersion and pho- 
tometric profiles agree with a constant mass-to-light ratio 
isotropic orbit model with M/Lb ~ 7, similar to what is ob- 
tained with anisotropic models; therefore it was concluded 
that dark matter is not needed within ~ 4.5Re. 

Romanowsky et al. (2003) found a decrease in the 



2 THE TOTAL MASS PROFILE OF NGC3379 

NGC 3379 is the dominant elliptical galaxy in the nearby 
Leo Group, and has an optical profile well described by a de 



^ This Is a ratio between spatial quantities, that is between 
the mass and luminosity enclosed within a sphere of radius 
r =Re (Saglia, private communication) . 
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Table 1. General and hot ISM projjcrtics of NGC.3379 



D 


Morph. 




Lb 


Re" 




band 




(Mpc) 


type 


(Bmag) 


(lOi"is,0) 


(", kpc) 


(km s (erg s 


(kcV) 


(kcV) 


10.6 


El 


10.18 


1.5 


37.5, 1.93 


221 4.6x10^7 


0.2-5 


0.52±0.30 



D is from Tonry at al. 2001; the morphological Type and Bj^ arc from RC3; derives from B"^ and the adopted distance. 

" Pclcticr ct al. (2000) 
central stellar velocity dispersion from Hallyday et al. (2001) 
Pukazawa et al. (2006) 



velocity dispersion profile of planetary nebulae extending 
out to 200". Within this radius, their isotropic modeling 
constrained M/Lb = 7.1 ± 0.6, a result consistent with 
no dark matter, given the possible Mt,/LB values (Ger- 
hard et al. 2001). Dekel et al. (2005) questioned this result 
with radially anisotropic modeling, and showed that a total 
M/ Lb ^ 15 at 5i?e is consistent with the data; this implies 
Mh/M» = 1.1 - 1.5 (for an assumed M^/Lb =6-7). 

Summarizing, the optical studies converge towards a 
very small dark matter content within Rc, of the order of 
MhlM-^ ^ 0.1, except for the Saglia et al. (1992)'s analysis 
that gives Mh/M* = 0.37 — 0.83. At larger radii, the dynam- 
ical modeling gives results consistent with no dark matter, 
but also with a global galactic value of Af^/M, — 1.1 — 1.5 
(Dekel et al. 2005) or even Mh/M, =0.7-5 (Sagha et al. 
1992). 

2.2 Results of X-ray studies and first reduction of 
the discrepancy 

The hot gas content of NGC3379 is very low, as known 
since ROSAT observations (e.g., Pellegrini 2005, Brown 
& Bregman 2001, Trinchieri et al. 2000). The most re- 
cent Chandra data have been used by F06 to derive the 
total mass within R^. From a parameterized gas density 
and temperature profile entered in eq. (1), F06 obtained 
M^t = 8.7 X lO"'Af0, with a quoted 'typical error' of 30% 
[therefore M^^^ = (8.7 ± 2.6) x lO^^A^o]. From within the 
same region, we derive a total M = 3.7 x 10^° Mq, by taking 
M/Lb at Re of Kronawitter ct al. (2000) (see Sect. 2.1) and^ 
0.42 X their quoted Lb- Therefore, we determine a discrep- 
ancy of a factor of ~ 2 between the X-ray and the optical 
estimate. 

F06 noted a larger discrepancy with the mass derived 

by Kronawitter et al. (2000), by a factor of ~ 7 (sec their 
Table 5); however, they likely underestimated the Kron- 
awitter et al.'s total mass within R^, for which they give 
just 1.3xlO^°M0. In fact, this mass cannot be smaller than 
{M/Lb)c X 0.42Lb = 3.3 x IO^Mq, where {M/Lb)c = 
5.6 is the central mass-to-light ratio of Kronawitter ct al. 
In any case, M,/Lb > 4 for the old stellar population 
of NGC3379 (Gerhard et al. 2001), which already implies 
M* > 2.5 X 10^° Af© within r = Re. 

All in all, the observational results that will be consid- 
ered in the following are that i) M^^t within Rc is a factor 
of ~ 2 larger than the total mass estimated by Kronawitter 

^ For the de Vaucouleurs (1948) law, the luminosity enclosed 
within a spherical radius r = iie is ~ 0.42 of the total one. 



et al. within the same radius; ii) when taking an average 
M^/Lb ~ 6 for the stellar population and the Lb value in 
Table 1, corresponds to Mh/Mt ~ 1.3 within R^. 

In the next Section, with the aid of hydrodynamical sim- 
ulations, we show how an overestimate of a factor of 1.5 — 2 
for the total mass of NGC3379 can be easily obtained when 
using eq.(l) instead of eq. (2), since the hot gas is not in hy- 
drostatic equilibrium, but instead it is outflowing over most 
of the galactic body. A summary of the hot gas properties de- 
rived by F06 is given in Table 1; even though resolved point 
sources have been subtracted before deriving these proper- 
ties, some residual emission from unresolved stellar sources 
is still contributing to the X-ray luminosity of the gas in 
Table 1. 



3 THE HYDRODYNAMICAL SIMULATIONS 

The time-dependent equations of hydrodynamics with 
source terms, and the numerical code used to solve them, 
are fully described in Ciotti et al. (1991). The model galaxy 
is initially devoid of gas, a situation produced by the galactic 
winds established by type II supernovae, early in the evo- 
lution of elliptical galaxies. Then the mass losses from the 
evolving stellar population start replenishing the interstellar 
space, and the resulting gas flow evolution is followed until 
an age of ~ 9 Gyrs is reached, that is the estimated age for 
NGC3379 (Terlevich & Forbes 2002 give an age of 9.01^^ 
Gyr). 

The X-ray luminosity in the energy band adopted by 
F06 has been calculated using the APEC code in XSPEC, 
for the density and temperature distribution of the gas given 
by the simulations (solar abundance heis been adopted) . The 
luminosity weighted emission temperature in the same en- 
ergy band has also been calculated for the models. 

3.1 The galaxy model 

A very good approximation of the de Vaucouleurs (1948) law 
is given by the Hernquist (1990) distribution, that has the 
advantage that all its dynamical properties can be expressed 
analytically. This distribution is 



where rc is the scale length and relates to the effective ra- 
dius of the de Vaucouleurs law as ~-Rc/1.8153. This law 
gives a description of the luminous mass distribution accu- 
rate enough for the treatment of hot gas flows, on the scales 
that are relevant for the problem (from few tens to several 
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thousands of parsecs). In the simulations we adopt a central 
grid spacing of 5 pc to allow for a better sampling of the 
inner regions. 

The radial density distribution of the dark haloes of 
ellipticals is not well constrained by observations; theoret- 
ical arguments favor a peaked profile (Ciotti & Pellegrini 
1992; Evans & CoUett 1997), and high resolution numeri- 
cal simulations of dissipationless collapse produce a density 
distribution oc near the center (Dubinsky & Carlberg 
1991, Navarro, Fronk, & White 1996). The model galaxy is 
therefore a superposition of two Hernquist density distribu- 
tions, one for the luminous matter p«(r) with rc = r* and 
M — M-t in eq. (4), and one for the dark matter ph(j) with 
M = Mh, ra = rh = fir,. 

In order to build a galaxy model where the masses of 
the two components and their distributions are consistent 
with observations and also produce the observed ctc, we solve 
and project the Jeans equations for a globally isotropic two- 
component Hernquist model. For this model we find that 
a very accurate fit of the projected, mass averaged stellar 
velocity dispersion within 7?e/8 is 



10' 



102 



103 



10^ 



GM, 



X 0.096 + 0.12 



Mh/M, 



with 



(3 = 2.8153 



Mh/M, 



1.8153, 



(5) 



(6) 



where {Mh/M,)e is the dark-to-luminous mass ratio within 
a radius r = Re. Thus for assigned Re, ctc, total Mh/M^ 
and {Mh /M,)c. wc derive Af, and (3 from eqs. (5) and (6). 
The typical stellar mass-to-light ratios resulting from this 
procedure, for {Mh/Mt)e = 0.1 — 1, are between 4 and 8. 
The model phase-space distribution function is everywhere 
positive (Ciotti 1996). 

3.2 Time evolving inputs 

The time evolving input ingredients of the numerical simula- 
tions are the rates of stellar mass loss from the aging stellar 
population, and the rate of SNIa heating. In the numerical 
code the exact mass return rate prescribed by the stellar evo- 
lution theory is used (see, e.g., Ciotti et al. 1991), updated 
to take into account the latest stellar population synthesis 
models (Maraston 2005) and more recent estimates for the 
mass ejection from stars as a function of their mass (Ren- 
zini & Ciotti 1993). The SNIa heating rate is parameterized 
as LsN{t) = EsN RsN{t) Lb, where Esn = 10®^ erg is the 
kinetic energy injected in the ISM by one SNIa, and the 
number of events as a function of time is _Rgjv(t)-£/s oc t'" , 
where the slope s describes the unknown decay rate (see 
also Pellegrini & Ciotti 1998). RsN{t) is normalized to give 
the SNIa's explosion rate in nearby E/SOs for t = 10 Gyrs, 
i.e., 0.18/i?5 SNu (where 1 SNu = 1 SNIa per 100 yrs per 
10^°Ls,o, hn = Ho/75), that has an associated uncertainty 
of ~ 30% (Cappellaro et al. 1999). Therefore we consider 
as a free parameter also the ratio i? between the adopted 
SNIa's rate and the best fit estimate of Cappellaro et al., 
with -& — 0.67 — 1.33. From a comparison between realistic 
modelings of the rate and its observation in high-« surveys 
currently available, a slope s ~ 1.0 — 1.2 has been favoured 
(Greggio 2005). 
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Figure 1. From top to bottom the radial profiles of the gas den- 
sity, temperature and velocity for representative models of inflow 
(solid line), partial wind (dashed line) and global wind (dotted 
line); see Sect. 3.3. 



Another heating source, though less important than 

SNIa's, is the thermalization of the stellar mass losses to the 
local "temperature" set by the stellar velocity dispersion. At 
each radius this is determined using the velocity dispersion 
profile obtained by solving the Jeans equation for the two- 
component Hernquist model in the globally isotropic case 
(Ciotti 1996). If anisotropy is present, this heating differs lo- 
cally from the value it has in the isotropic case. However, the 
difference is small for realistic anisotropy distributions, and 
the total heating over the galaxy is independent of orbital 
anisotropy, as a consequence of the virial theorem. In ad- 
dition, the heating due to thermalization of stellar motions 
is significantly smaller than SNIa heating, and therefore the 
effects of anisotropy on (r) are negligible [they are just 
indirect effects in eqs. (7)- (8) of the next Section, through 
the gas temperature and velocity terms]. 



3.3 Results 

In Sect. 3.3.1 below we adopt the underlying total mass dis- 
tribution on which most of the modeling of the optical data 
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converge (Sect. 2.1). In Sect. 3.3.2 we use as input the mass 
model found by F06 and also a "maximum" dark matter 
model. In both Sections we vary the SNIa's parameters i? 
and s to obtain X-ray gas properties as observed when an 
age of ~ 9 Gyrs is reached. When this is the case, wo apply 
to the radial distributions of the gas density and internal en- 
ergy, resulting from the numerical simulation, the standard 
procedure of eq. (1), in order to recover M^t(r) and com- 
pare it with the true underlying mass distribution M{r). In 
practice, wo calculate the "empirical' M^g^{r)/M{r) ratio 
from the output of the numerical simulations as 

Me^t(r) _ k dpT 
M{r) nrrippg dr ' 

This ratio is also given by eq. (3), where the V term can 

be explicitoly calculated from the momentum equation [see 
eq. (10) of CP04)]. In quasi-stationary conditions, as those 
reached by our models at the end of the simulations, one 
can neglect the du/dt term in the expression for V, and the 
full expression for eq. (3) becomes: 

M{r) ^ pg^ 2g dr ^ ' 

where p is the local mass injection rate due to stellar mass 
losses and u is the radial flow velocity. Eq. (8) allows for a 
quantitative evaluation of the separate contributions of the 
non hydrostatic terms to the M^gf{r)/M{r) ratio. 

By varying the input parameters, the resulting gas flow 
phases arc of a global wind, a partial wind (where the gas 
is inflowing only over a very small central region) and a 
global inflow; a sample of the radial profiles of gas velocity, 
temperature and density for representative models of these 
three phases at the present epoch is shown in Fig. 1 [see 
Pellegrini & Ciotti (1998) for a more detailed description of 
the gas flow evolution]. 



3.3.1 Optically derived total mass model 

In this experiment the galaxy model is based on the conclu- 
sions of Sect. 2.1: it has (Mft/M«)e = 0.13, in accordance 
with Kronawitter et al. (2000), and a global M^/M* = 1.5, 
consistent with the results of Dekel et al. (2005) and Saglia 
et al. (1992). The other parameters are set to the observed 
values given in Table 1. 

With this total mass profile, the hot ISM maintains in 
a global wind phase throughout its evolution, if the SNIa's 
rate is such that i9 = 1, regardless of the slope s. After a 
time equal to the age estimated for NGC3379, the hot ISM 
luminosity in the 0.2-5 keV band is less than derived by 
F06 by at least an order of magnitude. If the SNIa's rate 
is lowered, the hot gas fiow maintains in a wind phase over 
most of the galactic body, while it develops an inflow in a 
very small central region (<20 pc; see the partial wind case 
in Fig.l). In this way a higher luminosity is obtained at the 
present epoch, that is comparable to the observed one when 
1? is close to its lowest values {'& ~ 0.7). 

In Fig. 2 the derived M^j/Af is shown for a sample 
case of a partial wind with a hot gas luminosity of 1/(0.2 — 5 
keV)= 3.6 X 10^^ erg s~^ and an emission weighted temper- 
ature of 0.53 keV, both values close to those reported by 
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Figure 2. The radial trend of M^^^/M obtained from eq. (7), for 
the partial wind model described in Sect. 3.3.1 (solid line). Also 
shown are the terms on the right hand side of eq. (8): the radial 
derivative term (plus unity, dotted line), the term with the mass 

source (plus unity, dashed line), and M^^^/M according to cq. (8) 
(dot-dashed line). The dotted vertical line marks the position of 
Re- 

F06 (see Table 1)^ It is obtained for s = 1.1 and i? = 0.72, 
after a time of 8 Gyrs. M^t(r) overestimates M{r), by a 
factor of ^ 1.5 — 1.6 around Re- This factor is in the direc- 
tion of, and also enough for, the explanation for the resid- 
ual discrepancy between derived by F06 and the total 
mass derived with optical methods (Sect. 2.2). Note how the 
M^st /M profiles derived with eq. (7) and with eq. (8) main- 
tain remarkably close to each other over the whole numeri- 
cal grid. This means that the terms describing the deviation 
from hydrostatic equilibrium and entering eq. (8) are really 
accounting for the whole discrepancy measured with eq. (7) 
(and it is also an independent check of the reliability and 
accuracy of the numerical code) . 

3.3.2 X-ray derived and maximum dark mass model 

We study here a galaxy model with (M^/M,)e = 1.3, in 
accordance with the derived by F06 (Sect. 2.2); the 

global Mh/M, = 1.5 as in the previous Section 3.3.1, and 
the SNIa's parameters are varied. 

Again, for i9 = 1 global winds that are far less lumi- 
nous than the observed Lx are obtained, when the galaxy 
is 9 Gyrs old. Luminosity and temperature close to the ob- 
served values arc shown again when the flow is a partial 
wind, with approximately the same low SNIa's rates derived 

A contribution from unresolved sources to L(0.2-5 keV) in Ta- 
ble 1 (Sect. 2.2), and the uncertainty associated with this mea- 
surement, can account for the slightly lower L(0.2-5 keV) of the 
model. 
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Figure 3. The same as in Fig. 2, for a model with a larger fraction 
of dark matter content within Rg [(M^/Af,)e = 1.3], consistent 
with the X-ray analysis of F06 (see Sect. 3.3.2). Again i9 = 0.72 
and s = 1.1. 



Figure 4. The same as in Figs. 2 and 3, for the maximum dark 
mass model, where an inflow reverts into a PW at ~ 9 Gyrs (see 
Sect. 3.3.2). 



in the previous Section (i? ~ 0.7). This fact can be explained 

by considering that, although (Mh/M,)c increases from 0.13 
to 1.3, the total mass in the central region remains the same, 
and the absolute value of the binding energy per unit mass 
decreases by ~ 15% at the center; this is the outcome of im- 
posing fixed values for the central stellar velocity dispersion 
and Rc (sec cqs. 5 and 6). In addition, being Lb also fixed, 
the rate of mass and heat input are the same as in the set of 
models of Sect. 3.3.1. All in all, the net result is just a small 
reduction of the binding energy in the central region, where 
almost all of the emission is produced. 

Since the flow pattern is the same as in the previous 
Section, also the mass estimates suffer from the same 
problems. Fig. 3 shows such an estimate for a model with 
Lx = 2.5 X 10^^ erg s~^ and a temperature of 0.53 keV, 
after 8 Gyrs; the input SNIa's parameters are s = 1.1 and 
1? = 0.72. While similar to Fig. 2 in the general trends of 
the various quantites, the overestimate is even larger than 
in Fig. 2. Thus, even if (Mfe/M,)e were really as high as de- 
duced from assuming hydrostatic equilibrium, the derived 
M(^t would be wrong. This experiment also shows that the 
hot ISM is quite insensitive to the underlying mass distri- 
bution and partition between stellar and dark mass within 
Re, at least for these low dark matter contents. 

In the last experiment, we check whether reliable mass 
estimates can ever be obtained with the X-ray procedure 
of eq. (1), for this galaxy. This requires the hot gas to be 
roughly in hydrostatic equilibrium over most of the galaxy, 
as in a slow global inflow (e.g., Ciotti et al. 1991). In or- 
der to obtain an inflow over a large region, the dark matter 
content must be increased. Therefore, we assume the max- 
imum dark matter content allowed for by optical studies 
(Saglia et al. 1992), that is a global ratio of M^/M, = 5 and 



(Mh/M«)e = 0.83 (note that this latter value is also consis- 
tent with the estimate of F06, within the uncertainties; see 
Sect. 2.2). If 1? ~ 0.7, the galaxy indeed hosts a global inflow; 
however, after an age of 9 Gyrs, it is far more luminous than 
observed (Lx>10'"' erg s"^). For t? > 0.7, the flow pattern 
is again a global inflow during the first part of the evolution, 
but before reaching an age of 9 Gyrs, it slows down until it 
reverts its direction into an outflow, over most of the galaxy. 
The final phase is a partial wind, whose luminosity is much 
lower than in a model where the gas is still mostly inflowing. 
For example, if = 0.92 and s = 1.0, the transition epoch 
is at ~ 9 Gyrs, after which the luminosity decreases by 2 or- 
ders of magnitude in <1 Gyr. The observed L(0.2-5 keV) is 
reached between 10 and 10.5 Gyrs, with an average emission 
temperature of ~ 0.4 keV. Fig. 4 shows the mass profile that 
would be derived for this model: as for the previous cases, 
an overestimate is obtained at Re. This experiment shows 
that the X-ray procedure applied to NGC3379 is forced to 
produce mass overestimates, because the low observed X-ray 
emission requires that the gas is outflowing over most of the 
galaxy. 



4 CONCLUSIONS 

We have analyzed the discrepancy between optical and X- 
ray mass estimates of NGC3379, a well studied nearby ellip- 
tical for which optical investigations are consistent even with 
absence of dark matter, but for which recent X-ray studies 
derived a total mass much larger than derived with optical 
methods, within Rc- First we have reduced the discrepancy 
to a more plausible factor of ~ 2, by a reanalysis of the total 
mass within Re implied by previous optical studies. Then we 
have argued in favour of the hypothesis of lack of hydrostatic 
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equilibrium as the explanation for the residual inconsistency 
between the X-ray and optically derived total mass. 

With the aid of hydrodynamical simulations tailored 
specifically for NGC3379, we have tested whether hot gas 
properties similar to those observed can be recovered, and 
then what total mass is derived if the standard X-ray pro- 
cedure is applied to the flow resulting from the simulation. 
In a first experiment the total mass distribution underlying 
the gas flow is that derived by the modeling of the optical 
data, with {Mh/Mt)e = 0.13. For a SNIa's explosion rate 
consistent with current estimates, a hot gas flow with lu- 
minosity and temperature close to the Chandra values can 
be obtained. However, except for a very small central re- 
gion, the gas is largely outflowing, with a velocity increasing 
outwards. For such a flow, the X-ray method based on the 
hypothesis of hydrostatic equilibrium gives an overestimate 
of the total mass by a factor of ~ 1.6, within Re, and by >2 
at larger radii. 

In a second experiment, the underlying total mass dis- 
tribution is that derived by the modeling of the X-ray data 
within _Re, under the assumption of hydrostatic equilibrium; 
i.e., {Mh/Mt)e is increased to 1.3, while the global Mh/Mt 
is the same as before. Again, the observed hot gas properties 
are reproduced by a partial wind; therefore, the Eissumption 
of hydrostatic equilibrium and the use of the X-ray method 
produce a mass overestimate of a similar factor, showing the 
inconsistency of the whole procedure for this galaxy. This ex- 
periment also shows that the hot ISM is quite insensitive to 
the underlying distribution of dark mass within Re, at least 
for these low dark matter contents. 

Finally, wo have tried a "maximum" dark matter model, 
where the amount of dark mass is the maximum allowed 
for by optical studies [Mh/M, = 5 and (Afh/M.)c = 0.83, 
which is also consistent with the X-ray results] . The flow can 
now keep itself in a global inflow, the gas is close to hydro- 
static equilibrium over most of the galaxy, and the galaxy 
mass estimated with the X-ray method would be reliable, 
but the X-ray luminosity is orders of magnitude higher than 
observed. A luminosity as low as observed is obtained if the 
SNIa's rate is increased, so that the flow starts as a global 
inflow and becomes a partial wind in later times. Therefore 
the X-ray method again produces an overestimate of the 
total mass at Re. 

In conclusion, NGC3379 is a clear example where the 
departure from hydrostatic equilibrium of the hot ISM gives 
a consistent explanation for the discrepancy between the 
optical and X-ray-derived mass values. In particular, this 
departure (while also required to reproduce the observed low 
X-ray emission) accoxmts for M^^^ being an overestimate and 
also for the size of the overestimate. This analysis suggests 
that a pr(>vious understanding of the hot gas behavior is 
required before attempting dynamical estimates with hot 
gas-based studies. Here, for example, the very low hot gas 
content suggests by itself that most of the galaxy likely hosts 
a wind. 

Lastly, we can check whether the X-ray properties of 
NGC3379 can put any constraint in the controversy con- 
cerning the "dearth of dark matter in elliptical galaxies", 
for which NGC3379 has been suggested to provide an exam- 
ple (e.g., Romanowsky et al. 2003), later criticized by Dekel 
et al. (2005). We flnd that hot gas with an emission and aver- 
age temperature as observed by Chandra can be present for 



the whole range of dark mass amount and distribution sug- 
gested by optical works, just varying the SNIa's rate within 
its estimated uncertainty (i.e., it has to be not so high to sus- 
tain a global wind, for the lowest M^/M* values, and not 
so low that a global inflow lasts until the present epoch, for 
higher M^/Mf up to 5). Even higher dark matter contents 
than found by optical works can actually still correspond 
to partial winds with the required X-ray properties. This 
insensitiveness to the underlying mass is not totally new: al- 
though for much X-ray brighter galaxies, mass models con- 
strained by extended stellar dynamical data, and including 
different amounts of dark matter, were found to produce 
similar temperature and X-ray brightness profiles, within 
the uncertainties in the X-ray measurements (Bertin et al. 
1993). The situation may have improved with the Chandra 
and XMM— Newton data, at least for the X-ray brightest 
galaxies. 
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